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Abstract. From an analysis of the spectrum (4000 A to 
8800A) of HD 101584 it is found that most of the neu- 
tral and single ionized metallic lines arc in emission. The 
forbidden emission lines of [01] 6300A and 6363A and 
[CI] 8727A are detected, which indicate the presence of a 
very low excitation nebula. The Ha, Fell 6383A, Nal Di, 
D2 lines and the Call IR triplet lines show P-Cygni pro- 
files indicating a mass outflow. The Ha line shows many 
velocity components in the profile. The Fell 6383A also 
has almost the same line profile as the Ha line indicat- 
ing that they are formed in the same region. From the 
spectrum synthesis analysis we find the atmospheric pa- 
rameters to be T e //=8500K, log g=1.5, V tur b=13km s _1 
and [Fe/H]=0.0. From an analysis of the absorption lines 
the photospheric abundances of some of the elements are 
derived. Carbon and nitrogen are found to be overabun- 
dant. From the analysis of Fe emission lines we derived 



6100K±200 for the emission line region. 
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1. Introduction 

Humphreys and Ney (1974) found near-infrared excess in 
HD 101584 and suggested that it is a massive F-supergiant 
with an M-type binary companion star (Humphreys 
1976). However, HD 101584 (V=7.01, F0 Iape (Hof- 
fleit et al. 1983)) was found to be an IRAS source 
(IRAS 11385-5517) (Parthasarathy and Pottasch 1986). 
On the basis of its far-infrared colors, flux distribution and 
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detached cold circumstellar dust shell, Parthasarathy and 
Pottasch (1986) suggested that it is a low mass star in 
the post-Asymptotic Giant Branch (post-AGB) stage of 
evolution. 

CO molecular emission lines at millimeter wavelengths 
were detected by Trams et al. (1990). The complex struc- 
ture of the CO emission shows large Doppler velocities 
of 130 km s _1 with respect to the central velocity of the 
feature indicating a very high outflow velocity. Te Lintel 
Hekkert et al. (1992) reported the discovery of OH 1667 
MHz maser emission from the circumstellar envelope of 
HD 101584. The OH spectrum has a velocity range of 84 
km s _1 and shows two unusually broad emission features. 
Te Lintel Hekkert et al. (1992) found from the images ob- 
tained from the Australian Telescope, that the OH masers 
are located along the bipolar outflow. The post-AGB na- 
ture of HD 101584 is also suggested by the space velocity 
of the star derived from the central velocity of the CO 
and OH line emission. This velocity of V ra d — 50.3 ± 2.0 
km s _1 does not agree with the galactic rotation curve 
assuming it to be a luminous massive population I F su- 
pergiant. 

Bakker et al. (1996a) studied the low and high reso- 
lution ultraviolet spectra and the high resolution optical 
spectra of HD 101584. Based on the strength of Hel (see 
also Morrison and Zimba 1989) N II, C II lines and Geneva 
photometry, Bakker et al. (1996a) suggest that HD 101584 
is a B9 II star of T e// = 12000K ± 1000K and log g = 3.0. 
Bakker ct al. (1996b) also found small amplitude light and 
velocity variations and suggested that HD 101584 is a bi- 
nary with an orbital period of 218 days. 

The optical spectrum of HD 101584 is very complex 
and shows many lines in emission. In this paper we re- 
port an analysis of the high resolution optical spectrum of 
HD 101584. 
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2. Observations and analysis 

High resolution and high signal to noise ratio spectra 
of HD 101584 were obtained with the European South- 
ern Observatory (ESO) Coude Auxiliary Telescope (CAT) 
equipped with the Coude Echclle Spectrograph (CES) and 
a CCD as detector. The spectra cover the wavelength re- 
gions 5360-5400A, 6135-6185A, 6280-6320A, 6340-6385A, 
6540-6590A, 7090-7140A, 7420-7480A, 8305-8365A and 
8680-8740A. The spectral resolution ranged from 0.165A 
at 6150A to 0.210A at 8700A. We have also obtained 2.5A 
resolution spectra of HD 101584 from 3900A to 8600A 
with the lm telescope and UAGS spectrograph and a 
CCD as detector at the Vainu Bappu Observatory (VBO), 
Kavalur, India. In addition we obtained CCD spectra with 
the same telescope and Coude Echelle spectrograph, cov- 
ering the the wavelength region 4600A to 6600A with a 
resolution of 0.4A. All spectra mentioned above were used 
in this analysis. 

All the spectra were analyzed using IRAF software. 
The equivalent widths of lines were found by fitting a 
gaussian. For blended lines de-blending was done by fit- 
ting multiple gaussians. We carried out spectrum synthesis 
calculations using KURUCZ stellar models (1994). SYN- 
SPEC code (Hubeny ct al. 1985) was used for calculating 
the theoretical line profiles. The gf values were taken from 
Wiese et al. (1966), Wiese and Martin (1980), Hibbert 
et al.(1991), Parthasarathy et al. (1992) and Reddy et al. 
(1997 and references therein). For the analysis of forbidden 
lines we have used the IRAF software package NEBULAR 
under STSDAS. 

3. Description of the spectrum 

The remarkable characteristic of the optical spectrum of 
HD 101584 is the fact that different spectral regions re- 
semble different spectral types. The spectrum in the UV 
region is similar to that of a Lep which is an F-supergiant 
(Bakker 1994). The optical spectrum in the range 3600A- 
5400A is dominated by absorption lines. Most of them 
are due to neutral and single ionized lines of Ti, Cr and 
Fe. The Call H and K absorption lines are strong. The 
strength of the absorption lines are similar to that ob- 
served in an A2 supergiant. In the yellow and red spec- 
tral regions, most of the lines are in emission (Fig. 1). 
The emission lines show complex line profiles. The absorp- 
tion lines of NI, OI, CII and Sill are broad. The Paschen 
lines arc in absorption. Some of these absorption lines are 
blended with emission lines and many have asymmetric 
profiles. The OI lines at 6156A are blended with emission 
lines of Fel. The NI lines are strong and show asymmetric 
line profiles. The blue wing is shallow compared to the red 
wing. The CII lines at 6578A and 6582A are weak. The Na 
D lines, KI 7700 A (Fig. 2), the Call IR triplet lines (Fig. 
3), [OI], [CI] and Mgl 6318. 7A lines are found in emission. 
The OI triplet lines (Fig. 2) are very strong indicating an 
extended atmosphere and NLTE effects. 



3.1. P-Cygni profiles 

The Ha line has a very strong P-Cygni profile indicat- 
ing an outflow. The profile looks very complex. It shows 
at least 6 velocity components. The Fell line at 6383A 
is in emission and the profile is very similar to that of 
Ha (Fig. 4). Similar behaviour of the 6383A Fell line 
and Ha line is also noticed in the post-AGB F supergiant 
IRAS 10215-5916 (Garcia-Lario et al. 1994). The Ha and 
the Fell 6383A line show an outflow velocity of 100±10 
km s _1 . The H/3 line also shows a P-Cygni profile. It has a 
broad emission wing at the red end. This indicates that the 
line forming region is extended. The H/3, Nal Dl, D2 and 
the Call IR triplet lines (Fig. 3) show an outflow veloc- 
ity of 75±20 km s _1 . The velocity structure seen in these 
P-Cygni profiles could be due to emission from different 
shells formed during the episodic mass-loss events. 

3.2. Fel and Fell emission lines 

The presence of numerous emission lines of Fel and Fell 
makes it possible to derive the physical conditions of the 
line forming region. From the curve of growth analysis 
of the Fel and Fell emission lines (Viotti 1969), we have 
derived T ea;i =6300±1000K and 5550±1700K respectively 
(Fig. 5). The scatter found could be due to the fact that 
the lines are not optically thin. On the other hand, there 
are only few emission lines of Fell present in the spectra 
and thus the estimate from Fell might not be accurate. In 
order to determine whether the large scatter observed in 
Fig. 4 is reflecting optical thickness effects we have done 
self-absorption curve (SAC) analysis (Friedjung and Mu- 
ratorio 1987) for the Fel emission lines. 

SAC is a kind of curve of growth applied to emis- 
sion lines, but it has certain advantages as compared to 
the classical emission line curve of growth analysis. This 
method of analysis is valid also for optically thick lines. It 
deals with each transition separately, so that it is possible 
to get the population of different levels without assum- 
ing a Boltzmann distribution. In this curve, a function of 
the line flux emitted in the different transition of a given 
multiplet is taken in such a way that it is constant for a 
optically thin uniform medium. As the optical thickness 
increases the curve will move towards a straight line in- 
clined at -45°. The shape of the SAC in Fig. 6(a) shows the 
lines arc optically thick. The shape of the SAC is obtained 
by shifting all the multiplets with respect to a reference 
multiplet. Here we have taken multiplet 207 as reference. 
The X and Y shifts of each multiplet gives the relative 
population of the lower and upper level with respect to 
the reference multiplet. Fig. 6(b) and Fig. 6(c) shows the 
Y and X shifts versus the upper and lower excitation po- 
tential from which we derive the T e:r i=6100±200K. 
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3.3. Forbidden lines 

The forbidden emission lines at 5577A, 6300A and 6363A 
of neutral oxygen are present in the spectra. The forbidden 
line of neutral carbon at 8727A is also seen. The 6300A 
line is blended with Sell line and the 5577A line is very 
weak. We have calculated the I(6300)+I(6363)/I(5577) to 
be 13.3. From the flux ratio we can calculate T e (Oster- 
brock 1989). This flux ratio is not very accurate because 
of very the weak 5577A line and poor signal to noise spec- 
trum. For the flux ratio of 13.3 we derived a function de- 
pending on the electron density N e and temperature T e . 
Figure 7 shows the N e and T e contours for different values 
of flux ratio around 13.3. Since we do not see any other 
forbidden lines which are sensitive to the electron density, 
we could not fix both N e and T e uniquely. But assum- 
ing a temperature derived from the Fe emission lines, an 
electron density of 1 x 10 7 is obtained. For this value of 
electron density and temperature the C/O = 0.5± 0.2 has 
been obtained. 

4. Radial velocities 

There are very few absorption lines and most of these 
are affected by emission and or a shell component, there- 
fore we derived the average radial velocity from the well 
defined emission lines. The average radial velocity from 
the emission lines is found to be 50 ± 2km s _1 . Morrison 
and Zimba (1989) using 14 best absorption lines found 
the radial velocity to be 69± 1 km s _1 . From the equiv- 
alent widths of Fel absorption lines given by Rosenzweig 
et al. (1997) we find no correlation between log gf - x@ 
and heliocentric radial velocity (Fig. 8). However, Bakkcr 
et al. (1996a) found a correlations between log gf - \Q 
and heliocentric radial velocities of HD101584 in the UV. 
The discripancy could be due to the poor resolution of 
Rosenzweig et al. (1997) data compared to that of Bakkcr 
et al. (1996a). The large scatter seen in the radial veloci- 
ties could be due to pulsation. Similar velocity variations 
were noticed in other post-AGB supergiants (Garcfa-Lario 
et al. 1997, Hrivnak 1997). 

5. Atmospheric parameters and chemical 
composition 

The UV (IUE) low resolution spectrum of HD 101584 
matches well with that of an A6Ia star (HD 97534) (Fig. 
10) indicating a T eff of 8400K (Lang 1992). The pres- 
ence of CII lines at 6578 A and 6582A indicates a T e // 
> 8000K. For T e// < 8000K the CII lines would be very 
weak or absent. The Paschen lines also indicates a low 
gravity (Fig. 11). The luminosity class la also indicates 
a very low gravity. From the analysis of several nitrogen 
lines around 7440A and 8710A we derived the microturbu- 
lence velocity V tMr ;,=13km s _1 . We synthesised the spec- 
tral region from 4000A to 4700A (Fig. 9) with low gravity 
(log g =1.5 )models of Kurucz (1993) with temperatures 



8000K, 8500K and 9000K. The best fit was found for T eff 
= 8500K, log g = 1.5, V T =13km s" 1 and [Fe/H] = 0.0. 

The line at 5876A was identified as a Hel line by 
Bakker et al. (1996a) who also state that the lines at 5047A 
and 5045A as due to Hel and Nil respectively. However, 
we find that the 5047 and 5045 lines are in fact due to 
Fell. Except Hel 5876A, we have not found any other he- 
lium lines in the spectrum and nor have we found any Nil 
or Oil lines. In fact, Hibbert et al. (1991) indicate the 
presence of a CI line at 5876A. It is likely that the line at 
5876A may be due to CI instead of Hel. 

If we assume that the 5876A line is due to Hel then 
for a solar helium abundance and log g = 1.5, T e ff = 
9000K is found. Since we do not see any other helium 
lines, if 5876A line is due to helium , it is likely that it 
may be formed in the stellar wind or in the chromosphere 
of the star. On the basis of the presence of this helium line 
Bakker et al. (1996a) suggested that HD 101584 is a B9II 
star of T e ff 12000K. On the basis of the analysis of our 
spectra we have not found any evidence for such a high 
temperature. We have also analysed the equivalent widths 
of absorption lines in the spectrum of HD 101584 given by 
Bakker et al. (1996a). The final abundances of some of the 
elements are listed in Table 2. The abundances listed in 
Table 2 show that the star is overabundant in carbon and 
nitrogen. It appears that the material processed by the 
triple alpha C-N and O-N cycle has reached the surface. 

6. Discussion and Conclusions 

The optical spectrum of the post-AGB star HD 101584 
is rather complex. We find several emission lines and P- 
Cygni profiles indicating an ongoing mass-loss and the 
presence of a circumstellar gaseous envelope. From the 
analysis of the absorption lines we find the atmospheric 
parameters to be T e //=8500K , log g=1.5, V t =13km s _1 
and [Fc/H]=0.0. 

Carbon and Nitrogen are found to be overabundant 
indicating that material processed by triple alpha C-N and 
O-N cycles has reached the surface. Since our blue spectra 
arc of relatively low resolution and because of the presence 
of emission and shell components it is difficult to estimate 
reliable abundances of s-process elements. The OI line at 
6156A is blended with a weak Fel emission line. The OI 
triplet at 7777 A is very strong and affected by NLTE. In 
any case it appears that the oxygen abundance is nearly 
solar. A NLTE analysis of the high resolution OI 7777A 
triplet may yield a more reliable oxygen abundance. 

The nitrogen abundance is based on 6 lines in the 
7440A and 8710A region. We have not used the strong ni- 
trogen lines. Nitrogen seems to be clearly overabundant. 
The carbon abundance is based on two CII lines at 6578A 
and 6582A. There is a clear indication that carbon is over- 
abundant. The abundance of Mg, Ti, and Fe are nearly 
solar. The Ti abundance is based on 15 lines and the Fe 
abundance is based on 6 lines. Many of the other atomic 
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lines are affected by emission and shell components. In 
our opinion, the line at 5876A might be due to CI (Hib- 
bert et al. 1991) and not to Hel, as previously suggested 
by Bakker et al. (1996a). We have not found any other 
Hel, Nil or Oil lines. Our analysis shows that the T e ff is 
8500±500K. 

Bakker et al. (1996b) found small amplitude light and 
velocity variations and suggested that HD 101584 is a bi- 
nary with an orbital period of 218 days. The radial veloc- 
ity variations may be due to pulsation, macroturbulcncc 
motions or shock waves in the outer layers of the stellar 
atmosphere. Many post-AGB supergiants show small am- 
plitude light and velocity variations (Hirvnak 1997). These 
variations may not be interpreted as due to the presence 
of a binary companion. Long term monitoring of the ra- 
dial velocities is needed in order to understand the causes 
for these variations. 

The spectrum and the brightness of HD 101584 ap- 
pears to remain the same during last two or three decades. 
There is no evidence for significant variations in bright- 
ness similar to those observed in Luminous Blue Vari- 
ables (LBVs). The chemical composition and all the avail- 
able multiwavelength observational data collected during 
the last two decades by various observers indicates that 
HD 101584 is most likely a post-AGB star. 

The presence of several P-Cygni lines with significant 
outflow velocities, the OH maser and CO emission profiles 
(Tc Lintel Hekkert et al. 1992, Trams et al. 1990) and 
the IRAS infrared fluxes and colours (Parthasarathy and 
Pottasch 1986) indicates the possibility that HD 101584 
is a post-AGB star with a bipolar outflow with a dusty 
disk. Since HD 101584 shows a strong Ha emission line, 
high resolution imaging with the Hubble Space Telescope 
(HST) may reveal the bipolar nebula and the presence of 
a dusty disk similar to that observed in other post-AGB 
stars like IRAS 17150-3224 (Kwok et al. 1998) or IRAS 
17441-2411 (Su et al. 1998). 
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Fig. 1. High resolution spectra of HD 101584 obtained with the ESO CAT-CES 
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Fig. l(contd.). High resolution spectra of HD 101584 obtained with the ESO CAT-CES 
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Fig. l(contd.). High resolution spectra of HD 101584 obtained with the ESO CAT-CES 
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Fig. l(contd.). High resolution spectra of HD 101584 obtained with the ESO CAT-CES 
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Fig. 2. The spectrum in the upper panel shows several nitrogen lines and emission lines of Fe. The lower panel shows 
the KI 1699 A in emission and the strong absorption due to 01 triplet at 1111 A. 
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Fig. 3. Call IR triplet lines showing P-Cygni emission. This spectrum is of 2. 5 A resolution, obtained from VBO, 
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Fig. 4. P-Cygni profile of Ha and FeII(6383A) lines showing similar velocity structures 




Fig. 5. Curve of growth analysis of Fe emission lines. + represents the Fel lines and represents Fell lines. The 
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Fig. 6(a). 77ie p/ot shows the shape of the SAC. The + sign indicates multiplets 167, 168, 169,204-, 207, 208 °f Fe I having 
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Fig. 6 (b) . The fit shows the distribution of upper level population of different multiplets of Fel with respect to the 
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Fig. 6 (c) . The fit shows the distribution of lower level population of different multiplets of Fel with respect to the 
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Fig. 7. Plot of electron density N e and electron temperature T e . The dotted line is the contour for the observed ratio 
(13.3) of [OJI lines 5577A, 6300A and 6363A. Each contour in the plot is for a change in the flux ratio of O^j. 
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Fig. 9. IUE low resolution spectrum of HD 101584 is compared with that the A6Ia star HD 97534- The dotted line 
corresponds to the spectrum of HD 91534 while the solid line is HD 101584- 
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Fig. 10. Observed and synthetic spectra in the Paschen line region. 1-observed, 2-Teff=8000K,log g=1.0, 3- 
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Fig. 11. Synthesis spectra for different models are compared with the observed spectrum. The observed spectrum is of 
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Table 1 (contd.): List of emission lines detected in the high resolution spectrum 
of HP 101584 
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Table 1 (contd.): List of emission lines detected in the high resolution spectrum 



of HP 101584 
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Table 2: Chemical composition of HD101584. 



Element 


[element /H] 


Lines 




T e// =8500K 
log g=1.5 


T e// =9000K 
log g=1.5 




C 


1.0±0.1 


0.3±0.1 


C II 6578,6582 


N 


0.5±0.1 


0.6±0.1 


N I 7423,7442,8703, 
8711,8718,8728 


O 


0.1 


0.18 


O I 6158 


Mg 1 

Ti 1 

Fe 1 


0.2 

0.0±0.4 
-0.1±0.1 
0.2±0.4 


0.02 
0.4±0.4 
0.5±0.1 
0.3±0.4 


Mg II 4481.13 
Ti II 15 lines 
Fc I 3 lines 
Fe II 61ines 



